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In the adult hippocampus, neurogenesis — the process of generating mature 
granule cells from adult neural stem cells — occurs throughout the entire life- 
time. In order to investigate the involved regulatory mechanisms, knockout 
(KO) experiments, which modify the dynamic behaviour of this process, 
were conducted in the past. Evaluating these KOs is a non-trivial task 
owing to the complicated nature of the hippocampal neurogenic niche. In 
this study, we model neurogenesis as a multicompartmental system of ordin- 
ary differential equations based on experimental data. To analyse the results 
of KO experiments, we investigate how changes of cell properties, reflected 
by model parameters, influence the d3mamics of cell counts and of the exper- 
imentally observed counts of cells labelled by the cell division marker 
bromodeoxyuridine (BrdU). We find that changing cell proliferation rates or 
the fraction of self-renewal, reflecting the balance between symmetric and 
asymmetric cell divisions, may result in multiple time phases in the response 
of the system, such as an initial increase in cell counts followed by a decrease. 
Furthermore, these phases may be qualitatively different in cells at different 
differentiation stages and even between mitotically labelled ceUs and aU 
cells existing in the system. 



1. Introduction 

In the adult hippocampus, neurogenesis occurs in the subgranular zone of 
the dentate gyrus [1]. Here, individual stem cells are capable of generating 
astrocytes and neural progenitors. Recent data obtained from single cell level 
analysis demonstrate that stem cells perform four different types of events in 
order to produce progeny: symmetric divisions by dividing into two stem 
cells; two types of asymmetric divisions by either dividing into a stem cell 
and an astrocyte or a stem cell and a neural progenitor and astrogenic trans- 
formation; the direct conversion of a stem cell into an astrocyte [2]. Being 
born from stem cells, neural progenitors are capable of expanding their own 
pool by symmetric divisions and produce immature neurons called neuroblasts, 
which, in turn, mature to become neurons [2]. On a population level, it was 
shown that the number of stem cells, neural progenitors and immature neurons 
decreases during the ageing process [3,4] and alongside the number of newborn 
neurons depletes with time [5]. 

Over the past few years, knockout (KO) experiments have been used to study 
the neurogenic niche of the dentate gyrus, and the involved mechanisms govern- 
ing stem cells' fate choice [2,6-10]. Evaluating the results of such experiments is a 
non-trivial task owing to the multifactorial nature of the neurogenesis process. 
These complex dynamics severely limit intuitive interpretation of experimental 
data and call for tools such as mathematical modelling and analysis. The aim of 
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this study was to establish a basic model of adult hippocampal 
neurogenesis using previously published data. Because KO 
experiments had targeted stem cell compartments, a mathemat- 
ical model describing the d3mamics of cell counts for a given set 
of stem cell parameters provides a theoretical framework to 
identify the function of such KOs. 

Mathematical and computational models have been 
applied before to study adult neurogenesis. Ashbourn et ah 
[11] provide a system of partial differential equations to 
model the migration of immature neurons from the subventri- 
cular zone along the rostral migratory stream to the olfactory 
bulb and investigate parameters that lead to biologically 
plausible solutions. Aimone et al. [12] model the functional 
integration of new neurons to the hippocampus as an artificial 
neural network. To the authors' best knowledge, there exists 
no model addressing the cellular dynamics in the subgranular 
zone niche of the dentate gyrus. 

Our proposed model of the adult hippocampus is a neuro- 
genesis-adjusted modification of the model of haematopoiesis 
investigated by Marciniak-Czochra et al. [13] and Stiehl & 
Marciniak-Czochra [14]. Dynamics of hierarchical cell pro- 
duction systems, which maintain a continuous supply of 
differentiated functional cells to various parts of a living 
organism, have attracted the attention of biologists and 
mathematicians for many years in the context of blood cell 
production [15]. Besides common elements that can be found 
in all cell production systems, there are significant differences 
depending on the type of cells considered. To model the hierarch- 
ical structure of the system, we apply a system of ordinary 
differential equations (ODEs), each of which describes a discrete 
differentiation stage. In such models, the pace of commitment 
is dictated by successive divisions. However, in the case of 
neurogenesis, there are indications that stem cell differentiation 
also involves direct (continuous) transitions. Furthermore, 
neural stem cells are multipotent and generate, both, neuro- 
genic progenitors and astrocytes. We develop a new model 
accounting for these observations, as presented in §2. Another 
important application of modelling is in the choice of regulatory 
mechanisms. Because we aim to model short-term d3mamics 
of labelled cells, and there is no experimental evidence of 
feedback loops governing this process, we propose a linear 
model. This assumption stays in line with a parsimonious (reduc- 
tionist) approach to modelling, in which comprehensive models 
are better understood in view of simpler models. It allows closed- 
form solutions to be obtained for the mathematical analysis of 
derivatives with respect to stem cell parameters. 

Our study is organized as follows: in §2, we state an ODE 
model of adult hippocampal neurogenesis based on the 
experimental observations reviewed in the first paragraph of 
this introduction. Moreover, we introduce parameters that 
model the dynamics of neural stem and progenitor cells, 
namely the fraction of self-renewal, the proliferation rate 
and the division probability. In §3, we infer relations among 
these model parameters by deriving parameter conditions 
that account for the age-related decline in stem cell and pro- 
genitor counts as demonstrated by experimental data. 
Section 4 provides a mathematical analysis of the effects of 
a KO experiment. Because a stem-cell-targeting inducible 
KO spontaneously changes the dynamics of its target, we 
model such a KO by analysing the effects of alterations (calcu- 
lating partial derivatives) with respect to the stem cell 
parameters proliferation rate, fraction of self-renewal and div- 
ision probability on cell counts and on the number of 



bromodeoxyuridine (BrdU) incorporating cells. Section 5 con- 
tains parameter estimations and numerical investigations that 
could not be treated analytically and, in §6, we summarize 
and discuss our findings. 

Basic notation: we occasionally write x{t; p) to emphasize the 
dependence of the solution x{t) of a differential equation on a 
parameter p and sgn(a) denotes the sign of a real number a. 



2. Derivation of a multicompartmental model 

Based on the experimental data outlined in §1, we assume a 
model of adult neurogenesis with five cellular compartments: 
stem cells (ci), (neural) progenitors (C2), neuroblasts (C3), 
mature neurons (q) and astrocytes (C5): 

^ = {2aiei-l)pici{t\ 
dc? 

= 0i2(l - ai)KpiCi(t) + ((2fl2^2 -1)P2- d2)C2{t), 

dc^ 

— ^ = ((1 + f?2) - 202^2) plClit) - p3C3it) - d3C3{t), 

= P3C3(t) - diCiit), 
^ = (#12(1 - fli)(l - k) + 1 - 0i)piCi(f) - rf5C5(0- 

(2.1) 

The model equations describe the following events: stem cells 
can either divide with probability 9i or perform an astrocy- 
tic transformation with probability \ — 9i. The rate at which 
stem cells undergo these events is given by the parameter 
pi, called hereafter as the proliferation rate. The balance 
between symmetric self-renewal and asymmetric divisions 
is reflected by the fraction of self-renewal a^ > 1/2, which is 
the probability that a daughter cell has the same fate as the 
mother cell. The relation 

ai = l(2ai-l) + ^(l-(2ai-l)), 

shows that 2ai - 1 is the corresponding probability of a sym- 
metric division and 1 - (2ai - 1) = 2(1 - ai) the probability 
of an asymmetric division, because the fraction of daughter 
cells that continue as stem cells is 1 in a symmetric and 1/2 
in an asymmetric division. It follows that the expected net 
change of the number of stem cells after one stem cell event 
(division or transformation) is given by 

Oiilai - 1) + (1 - ^i)( - 1) = laiOi - 1. 

Asymmetric cell divisions may lead to two types of differen- 
tiated cells. The non-stem daughter cell is assumed to be 
either a neural progenitor with probability k or an astrocyte 
with probability 1 — k (see figure 1 for the diagram showing 
possible scenarios followed by a stem cell). 

For the proliferative capacity of progenitors, we again 
assume two possible modes of generating progeny: division, 
which occurs with probability O2 or direct transformation to a 
neuroblast with probability 1 — 62. Analogous to stem cells, 
progenitors have a proliferation rate p2 and a fraction of self- 
renewal ^2- Neuroblasts mature by transforming into a neuron 
with rate p^. Furthermore, we assume that aU cell t3rpes 
except stem cells are subjected to decay via apoptosis, modelled 
by the parameters di corresponding to cellular compartment i. 
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Figure 1. Proliferation diagram of a stem cell. Red nodes indicate events with stochastic outcome (e.g. division or transformation; symmetric or asymmetric division), 
blue nodes describe the outcome of particular events using chemical reaction notation (S, stem cell, P, neural progenitor. A, astrocyte). 0^ denotes the probability of 
stem cell division, pi denotes the proliferation rate; ^7^ reflects the probability that a daughter cell has the same fate as its parent cell (self-renewal takes place) and 
K is the probability that a neural progenitor is produced in an asymmetric division rather than an astrocyte. 
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The independent time variable t is used in two contexts. In 
§3, we analyse age-related properties of the neurogenesis 
system and use time for the adult age of the animal, i.e. the 
time point, t = 0, refers to the beginning of adult age, and the 
initial data consist of the number of cells present at f = 0 for 
each compartment. In §§4 and 5, we investigate the effect of 
altered stem cell dynamics on cell counts and on the number 
of BrdU-labelled cells in the framework of an inducible KO 
experiment, which imposes altered d3mamics upon adminis- 
tration of a certain chemical. In this context, time is used as the 
time after the KO occurred, i.e. t = 0 refers to the start at which 
the neurogenesis system operates under the altered dynamics. 



3. Decline of stem cell and progenitor counts 

Because it was observed that the number of stem and pro- 
genitor cells declines with age [3], we first derive parameter 
conditions that account for this effect. The proof of the 
subsequent lemma is given in appendix A 1. 

Lemma 3.1. The solutions Ci{t) and C2{t) of (2.1) are monotoni- 
cally decreasing for all t>0 if and only if 

1 



(^lOi < - and 



ci(0) ^ d2-{2a202 -l)p2 
C2(0) K0i2(l-ai)pi 



Furthermore 



lim^ f > 0 {2ai0i -l)pi> (2^202 - 1)^2 - ^^2, 
t^oo C2(t) \ = 0 otherwise. 



Biological interpretation: Lemma 3.1 states that the deple- 
tion of the stem cell pool takes place if and only if symmetric 
stem cell divisions, accompanied by a gain of stem cells, are 
less likely than astrocytic transformations with the resulting 
loss of the stem cell. The second part states that the ratio 
of the number of stem cells to the number of progenitors 
converges to zero, if the net depletion rate of stem cells is 
higher than the one of progenitors. Otherwise, it goes to a posi- 
tive value. Furthermore, the positive steady state is attained 
monotonically, either increasing or decreasing. Interestin- 
gly, both behaviours have been observed experimentally. 
In reference [3], it was reported that the ratio of the number 
of stem cells to the number of progenitors is monotonically 



decreasing, whereas Jinno [16] reports an increasing pro- 
gression. The discrepancy in both observations might thus 
come from different labelling paradigms and observations of 
different subpopulations. 

In [3, electronic supplementary material, table S2], a time 
series for the age-related decline of the stem cell and progeni- 
tor count was obtained. Fitting these data to the solution of 
(2.1) indicates that the parameters of our model satisfy the 
relations (2^2 ^2 ~ l)/^2 ~ ^^2 ~ (2^1 ^1 - l)pi < 0 and aiOi < 
1/2 (see §5.1). Hence, these data are consistent with the scen- 
ario in which the net depletion rate of progenitors is higher 
than the net depletion rate of stem cells. To summarize, we 
further make the following assumption in our subsequent 
mathematical analysis in order to remain consistent with 
the experimental data of reference [3]: 

Assumption 3.2. The parameters of model (2.1) have the 
properties 

(2fl2^2 - '^)P2 -d2- (2ai0i - l)pi < 0, 

and 

1 

diOi < -. 



4. Mathematical analysis of altered stem cell 
parameters 

4.1. Preliminaries 

4.1.1. The inducible knockout experiment 

A gene-KO is a procedure that eliminates a certain gene (a 
DNA sequence encoding a protein) from an organism's 
DNA. Thus, the corresponding protein is not synthesized. 
In order to study the protein-driven regulatory mechanisms 
involved in the dynamics of adult neural stem cells, such 
KO experiments were conducted in the past. A particular 
version of KO experiments is the inducible KO: cells with 
a pre-marked gene-sequence react to the administration of a 
chemical that is injected in the animal, with the activation of 
a cutting enzyme that excises the marked sequence. Because 
it is possible that this cutting enzyme is only present in cells 
expressing a certain other protein, for instance, the adult 
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Figure 2. Graphical representation of the analysed knockout (KO) scenarios, (a) Cell counts C/(f) are evaluated at time t after the KO. (b) The number of BrdU- 
labelled cells is evaluated at time r after the end of the labelling period (r = 0) and BrdU was given at time f* after KO. 



stem cell characteristic protein nestin, one can selectively 
knock out the gene of interest in stem cells. 

If the KO of a certain gene in the stem cell compartment 
results in a difference between KO and wild-type (non-KO) 
animals regarding the number of counted cells, the question 
arises as to which stem cell parameter was affected by the 
KO and caused the observed difference. To treat this ques- 
tion in a general way, we examine the effects of alterations 
of the stem cell parameters Ui (fraction of self-renewal), 6i 
(division probability) and (proliferation rate) on the 
number of existing cells and the number of cells labelled 
by BrdU, a chemical that is incorporated in the cells' DNA 
after cell division has taken place and during the stage of 
DNA synthesis. 

4.1.2. Modelling two experimental scenarios 

We consider two scenarios related to KO experiments, for 
which we analyse the effect of altered stem cell parameters. 
Scenario (i) — (figure la): starting from a time point zero, 
which corresponds to the fixed age of the studied animal at 
which the KO is conducted, the number of cells of compart- 
ment i is analysed at t time units after time point zero. 
Scenario (ii) — (figure 2h): at time after the initial time 
point of the KO, BrdU is administered and is present in the 
organism for a duration 5, thus labelling DNA-synthesizing 
cells during that period. At r time units after the labelling 
has ended, the number of BrdU-labelled cells (//(f*, r)) 
is examined. 

To evaluate for both scenarios the effects of a change of a 
parameter p G {ai, pi, Oi} from a value p to a value p -\- ^p 
{Ap > 0), we analyse the sign of the derivative dpCi{t; p), 
respectively, dpli{t^, r; p) with respect to the parameter p. 
Thus, we assume that Ap is so small, that 

sgn{ci{t; p + Ap)- Ci{t) p)) - s^n{dpCi{t; p)), 

and 

sgn(/,(f*, T,p + Ap) - kit'', r; p)) = sgn(ap/,(f*, r; p)). 

To model scenario (i), we use our model (2.1) together 
with initial data C/(0) corresponding to the number of cells 
of compartment z, which are present at the time point of 
the KO. For scenario (ii), the initial data are //(f*, 0), the 
number of cells that have incorporated BrdU at the end of 
the labelling period, where BrdU was given at time point 

and the independent variable in this scenario is r, the 
time that passed since the end of the labelling period. 

4.1.3. Initial data for BrdU-labelled cells 

It is known from the theory of branching processes that 
in a model of proliferation in which particles (cells) have 
exponentially distributed lifetimes, mean counts of particles 
(cells) follow a system of ODEs [17, ch. 4]. Conversely, if the 



population of cells is described by a system of ODEs, then 
the interpretation is that the cells have exponentially distrib- 
uted lifetimes. This is a simple and widely used model (see 
relevant discussion of the cell proliferation models in 
Kimmel & Axelrod [17]), despite the fact that it has been 
known that cell lifetime distributions are not exponentially 
distributed [18]. Thus, we use this relationship between 
ODEs and the exponential distribution in order to derive 
equations for the number of cells belonging to type / and 
have incorporated BrdU (//(t*, 0)). 

Recall that BrdU is a chemical that is incorporated in cells 
after they performed a division and are in the stage of DNA 
synthesis. Thus, in order to be labelled by BrdU, a cell must 
be in S phase during the time interval of length 8 in which 
BrdU is present in the animal. Assuming that the fraction 
of dividing cells equals the fraction of DNA-synthesizing 
cells during any time interval of fixed length, it follows that 
the number of cells that have incorporated BrdU at the end 
of the labelling period, starting with labelling at time f*, 
is given by 



b 



hit*, 0) : 

hit*, 0) : 

hit*, 0) : 

kit*, 0) : 

hit*, 0) : 



:2aiei{l-e-P'')ci{t*), 
:/<f?i2(l-fli)(l-e-''i«)ci(f*) 

+ 62i2a2-m-e-P'')c2it*), 
: e22il - fl2)(l - e-i^')c2it*), 

0, 

:(l-K)0i2(l-fli)(l-e-'"«)ci(f*), 



(4.1) 



where Ci{t^) is the number of cells of compartment z, present 
at time of BrdU injection. More specifically, a random 
variable X, exponentially distributed with parameter A, has 
the property 



Pr({X G [f, f + 8]]\{X >t]) = l-e- 



(4.2) 



Thus, the fraction of stem cells performing a transformation 
or division during marker exposure is 1 — e~P^^, and a frac- 
tion 01 of them divides with each division contributing, on 
average, 2ai stem cells. Analogous considerations lead to 
the equations for all other cell compartments. Note that neur- 
ons are assumed to be the result of a transformation from 
neuroblasts rather than a division. Thus, there are no BrdU- 
labelled neurons right after the BrdU-labelling period 
has ended. 



4.2. Effects of altered stem cell parameters 

Based on the considerations of §4.1.3, we model the scenarios 
(i) and (ii) of §4.1.2 using the system of ODEs that follow. The 
equations for the quantities //(f*, r), describing the number 
of BrdU-labelled cells of type i at r time units after BrdU 
was started at time t = t^, have been derived based on the 
assumption that labelled cells of type i follow the same 



dynamics as their corresponding compartment Cj. Thus, the 
initial data for these labelled cells, //(f*, 0), depend on 
C/(f*), the number of cells present at the time point t = 
of BrdU injection. 



dci 
Id 
dc2 
'dt 
dcs 
dt 

dC4 

'dt 

dC5 

dt 

Ci{0) = ni, 



{t) = {2aiei-l)pici{tl 

(t) = ^i2(l -ai)KpiCi(0 + ((2a2^2 - I)p2-d2)c2{t\ 
(t) = ((1 + ^2) - 2a2 02)p2C2{t) - dscsit), 

{t) = p3C3(t)-d4C4(t), 

(t) = (0i2(l -ai)(l - k) + 1 - 0i)pici(t) - d^c^{t\ 



f^(f*, r) = 01 2(1 -ai)KpMt*, t) + ((2^262 - l)p2 -rf2)/2(f*, r) 

QT 



(t*, t) = ((1 + »2)-2a202)p2kit*, r)-dMi*, r), 



{t*,r) = p3h{t*,r)-dik{t*,r), 



{t*, t) = (f?i2(l - k) + 1 - (h)pMt*, T)-dMt*, r), 



dr 
dr 

/i(f*,0) = 2fliei(l-e-'"Vi(f*), 

hit*, 0) = /c0i2(l -e-'"«)ci (t*)+2fl2 02(1 -e-''^«)c2(t*), 
/3(f 0) = 022(1 -fl2)(l -e-f'')c2{t*), 
h{t*,Q) = Q, 

/5(f*,0) = (l-/c)ei2(l-fli)(l-e-P'Vi(f*), 

(4.3) 

where Ui is the number of cells of compartment / at the fixed age 
of the animal where the KO is induced. We mathematically 
analyse the derivatives of Ci{t), C2{t) and li{t, 0), hit^O) and 
IsitfO) with respect to the parameter p G {ai, 6i, pi], thus 
evaluating the implications of altered stem cell parameters on 
the number of existing stem cells and progenitors and on the 
number of BrdU-labelled stem cells, progenitors and astrocytes 
at the end of the labelling period. The parameter-derivatives of 
the number of neurons and BrdU-labelled neurons (c4(f) and 
li{t,T)) cannot be analysed in full generality for arbitrary r 
and are thus investigated numerically in §5. For notational 
convenience, we define 

/f(f*) := li{t*, 0), 

and denote this quantity as the number of BrdU incorporating 
cells of compartment i. 

By a straightforward calculation (see appendix A 2), one 
obtains the identities 



daMi'rai) = 2eipitci{t) 
Op^cift pi) = (laiOi - l)tci{t) 
de^ciit; Oi) = 2aipitci(tl 



and 



dpCi{t;p) 



ci{t\ 



(4.4) 



(4.5) 



In order to evaluate (4.5), we need a technical lemma. 

Lemma 4.1. For t >0, a < 0 and /3, y G [R, consider the 
function 

_ {e-' -l){ap+y)-ayt 

It holds /(O) = 0, f{t) G 0{t) and f has the following 
properties: 

(PI) J/ /3 > 0 and 7 > 0, 

/ is monotonically increasing: 

(P2) If /3 > 0 and y < 0, there exists a unique tQ> such 
that for all t>0 

sgn(f(0) = -sgn(^ - to). 

(P3) If j8 < 0 and y > 0, there exists a unique to> 0 such 
that for all t>0 

sgn(f(0) = sgn(f-fo). 

Proof Follows from evaluating f{t) = {{ap + y)e'^^ — y)/a for 
every single case. 

Furthermore, we introduce the notion of the sign-sequence 
of a real-valued function: 

Definition 4.2. (sign-sequence). Let /: [0, oo) ^ [R be a func- 
tion. The sign-sequence a{f) is defined as the sequence of 
distinct signs of f{t) that are encountered by traversing the 
domain of / from zero to infinity. For instance, 
(7((f-l)(f-2)) = ( + ,0, -,0, +). 



4.2.1. Altered fraction of self-renewal 

Here, we describe the effect of increasing the fraction of 
self-renewal of stem cells. Recall that this fraction ai is 
defined as the probability that a daughter cell, which 
resulted from a stem cell division, becomes a stem cell 
itself. Thus, increasing the fraction of self-renewal increases 
the proportion of symmetric stem cell divisions that give 
rise to two stem cells at the expense of asymmetric stem 
cell divisions. 

Lemma 4.3. The solution of (4.3) satisfies 



o-(8aiCi(f;ai)) = (0, + ), 
cT{^aA{t;a^)) = { + \ 
o-(8aiC2(^;ai)) = (0, - , 0, + ; 

a(daMt;ai)) = {-,0, +). 



(4.6) 



for ; e [2,5] and p e jfli, Oi, pi), where a, p and y depend on 
both, / and p (cf. table 2). 



Proof. From (4.3) and (4.4), it follows that 
daiCi(t; ai) = lOipitciit), 
d„J°it; fli) = 2^1(1 - e-'"«)(l + 2fli0i pit^t), 
daM «i) = KOiil - €-'"«)( - 2 + 2(1 - fli)20i piOci(f) 

daM fli) = (1 - K)e^{l - e-'"«)( - 2 + 2(1 - fli)2ei pif)ci(f). J 

(4.7) 

Thus, da^Ci{t; ai) and da^Ci{0; ai) are positive for all positive t 
and daiCi{0; ai) = 0. Consider now the quantity da^C2{t;ai). 
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Figure 3. Simulated responses to an infinitesimal increase in stem cells fraction of self-renewal of the number of stem cells (c^), progenitors (cii, BrdU 
incorporating stem cells (/f) and BrdU incorporating progenitors (/2), respectively, at time f after the increase. The parameter set of §5.2 was used. 



Assumption 3.2 together with table 2 implies that the first 
factor of da^C2{t;ai) in (4.5) has the property (P3) in lemma 
4.1. Consequently, da^C2{t; ai) is negative on (0, to) and posi- 
tive on {to, 00 ) for some positive on to. From (4.7), we 
deduce that the same is true for daJ2i^)Cii) and daj^{t;ai), 
because the term - 2 + 2(1 - ai)20ipit is negative for t = 0 
and positive for sufficiently large t. ■ 

Biological interpretation. Lemma 4.3 shows that increasing 
stem cells' self -renewal increases the stem cell count and the 
number of BrdU incorporating stem cells at any time 
point after the increase was performed. Conversely, this 
increased self-renewal initially decreases the progenitor 
count and the number of BrdU incorporating progenitors 
and astrocytes. But this decrease is reversed and turns into 
an increase after the initial period. Furthermore, the effect 
of altered self-renewal is instantaneous on the number of 
BrdU-labelled cells in the sense that the corresponding 
parameter-derivative is non-zero at time zero. Figure 3 depicts 
a simulation of the time-dependent responses of an increased 
fraction of self-renewal. This simulation is consistent with 
lemma 4.3. 

The two-phase progression on progenitors and astrocytes 
can be explained intuitively as follows: the increased number 
of symmetric stem cell divisions at the expense of asymmetric 
divisions reduce the proportion of events at which progeni- 
tors are bom. Thus, a decreased number of progenitors is 
observed initially. At the same time, the increased number 
of symmetric stem cell divisions, which result in an enlarged 
stem cell pool, benefits the progenitor count in the long run: 
although a reduced fraction of stem cells generates progeni- 
tors via asymmetric divisions, the increased number of 
stem cells dominates this effect, meaning that the total 
number of asymmetric stem cell divisions is elevated. The 
immediate effect on BrdU incorporating progenitors and 
astrocytes can be explained by the observation that changing 
a parameter that affects division instantaneously changes the 
output of the division and that the mitotic marker BrdU 
exactly labels this output. Thus, the number of cells labelled 
by BrdU is faster influenced by a parameter change than 
the actual cell count. 



4.2.2. Altered proliferation rate 

The proliferation rate pi is the rate at which stem cells 
undergo division or transformation events. Increasing this 
rate shortens the waiting time between successive events of 
a given stem cell. 



Lemma 4.4. The solution of (4.3) satisfies 
cT(d,,Ci(t;pi)) = (0, -), 

(7(ap,/?(f;pi)) = ( + ,0, -), 
(T{dp,C2{t; pi)) = (0, + , 0, - ), K (4.8) 

a{d,f^{t;pr)) = {^,0, - ), 

a{d,f,{t;pr)) = {^,0, - ). 



Proof Equations (4.3) and (4.4) imply 

dp,ci(t; pi) = (2ai0i-l)tci(t\ 

dpf,{t; pi) = 2aiei{8e-P^' + (1 - e-P^'XlaiOi - l)t)ci{t\ 



ap,/^(f; pi) = K^i2(l-ai)(5e-Pi^+(l-e-Pi^)(2ai^i - l)t)ci{t) 

^2a2e2{l-e-P'')dpMi) Pi\ 
dpf,{t; pi) = (1 - k)0i2(1 - ai){8e-P'' 

+ (l-e-Pi^)(2ai0i-l)f)ci(O. 



(4.9) 



Because of the assumption 3.2, 2ai6i — 1 < 0. Thus, 
QpiCi{t; pi) is negative on (0, oo) and dpj^{t; pi) is positive 
on (0, to) and negative on {to, oo) for some positive to. The 
same is true for dp^l^{t; pi), but with different to. Consider 
now dp^C2{t; pi). The first factor of dp-^C2{t; pi) in (4.5) has 
the property (P2) in lemma 4.1. Hence, dp^C2{t;pi) and 
dpJ2{i; pi) show the same qualitative progression as stated 
for dpj^{t; pi). Analogous to the considerations of an altered 
fraction of self-renewal, the BrdU incorporating quantities 
dpJiit; pi),dpJ2{i; pi) and dpj^{t;pi) are positive at time 
zero, whereas dp^Ci{0; pi) = dp^C2{0; pi) = 0. ■ 

Biological interpretation. Lemma 4.4 shows that an increased 
proliferation rate of stem cells decreases the stem cell count, 
whereas the number of BrdU incorporating stem cells initially 
increases and later on decreases. This effect can be explained 
by the observation that an increased proliferation rate causes 
more stem cell divisions over any given time interval, resulting 
initially in more BrdU incorporating stem cells. Furthermore, 
the increased proliferation rate also causes a higher rate of astro- 
cytic transformations, the events responsible for the depletion of 
the stem cell pool. As time progresses, the increased decay rate 
of stem cells compensates for the higher proportion of BrdU 
incorporating cells and results thus in a net decrease of labelled 
stem cells. See figure 4 for a corresponding simulation, which is 
consistent with lemma 4.4. It is not surprising that the number 




Figure 4. Simulated responses to an infinitesimal increase in stem cells proliferation rate pi of the number of stem cells (Ci), progenitors [Cji, BrdU incorporating 
stem cells (/f) and BrdU incorporating progenitors respectively, at time f after the increase. The parameter set of §5.2 was used. 
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Figure 5. Simulated responses to an infinitesimal increase in stem cells division probability 0^ of the number of stem cells (ci), progenitors (q), BrdU incorporating 
stem cells (/f) and BrdU incorporating progenitors (If}, respectively, at time f after the increase. The parameter set of §5.2 was used. 



of progenitors and BrdU-labelled progenitors and astrocytes 
display the same qualitative trend as labelled stem cells, because 
these quantities also depend on stem cell divisions. 

4.2.3. Altered division probability 

The division probability Oi of a stem cell is the probability 
that the next event a stem cell undergoes is a division 
rather than a transformation. Consequently, increasing the 
division probability causes more division and fewer trans- 
formation events. 

Lemma 4.5. The solution of (4.3) satisfies 







= (0, + ), ^ 


a{de 




= ( + ), 




C2{t;ei)) 


= (0, + ), > 


<T(de 




= ( + ), 






= ( + )• , 



Proof It holds 
de,ci{t;0i) = 2aipitci{t), 
deAt; Oi) = 2ai{l - e-P''){l + Oilui pit)ci{t\ 
deA^-^ Oi) = k2(1 - ai)(l - e-P'')(l + Oilai pit)ci(t) 

deAt; Oi) = {1- k)2(1 - ai)(l - e-P'')(l + Oilai pit)ci(t). ^ 

(4.11) 

Thus, deiCi(t; 61), deJi{t; Oi) and deAt; Oi) are positive for 
positive t. Furthermore, the first factor of d0^C2{t; 61) in (4.5) 
has the property (PI) in lemma 4.1. Hence, d0^C2{t;6i) and 
deJ2i^) Oi) are also positive for positive t. ■ 

Biological interpretation. From lemma 4.5, we conclude that 
increasing the stem cell division probability causes an 



increase in cell counts and in the number of BrdU incorporat- 
ing cells for all considered compartments. There is no two- 
phase progression displaying increased or decreased cell 
numbers in the distinct phases after the change in stem cell 
dynamics. The effect on the progenitor count is qualitatively 
the same as on the stem cell count; the same holds true for 
BrdU incorporating stem cells and progenitors. Figure 5 illus- 
trates a corresponding simulation of the discussed quantities. 

4.3. Final remarks 

The above considerations show that all derivatives of the 
quantities Ci, C2, I2 and ^ with respect to the stem cell par- 
ameters fraction of self-renewal, proliferation rate and 
division probability are products of the exponentially 
decreasing function Ci{t) and a factor that can be bounded 
by an affine linear function of t. Thus, the effects of an altered 
stem cell parameter on the number of existing cells and the 
number of BrdU incorporating cells weakens with time 
owing to the decline of the stem cell compartment. The 
effect of altered stem cell parameters on the number of exist- 
ing astrocytes, dpC5{t;p), could not be analysed owing to no 
available data on the sign of af3 -\- y with a, ^ and 7 as 
stated in table 2. 

The sign sequences of the parameter-derivatives of the 
five quantities stem cell count, progenitor count and BrdU 
incorporating stem cells, progenitors and astrocytes with 
respect to the three considered parameters ai, pi and 61 are 
summarized in table 1. 



5. Numerical investigations 
5.1. Parameter estimations 

Our neurogenesis model (2.1) involves 18 parameters, 
including five parameters for the initial data of each compart- 
ment. We first analyse the parameter region that is consistent 
with the data presented in [3, electronic supplementary 



Table 1. Time-dependent responses of stem cells (c^), BrdU incorporating stem cells (/f), progenitors {cj), BrdU incorporating progenitors (Ij) and BrdU 
incorporating astrocytes (/5), respectively, to an infinitesimal increase of the respective parameter p. 
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Figure 6. Fitting model (2.1) to the number of stem cells and the ratio of the number of stem cells to the number of progenitors in [3, electronic supplementary 
material, table S2] results in a poor agreement with the number of progenitors. 
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Figure 7. Fitting model (2.1) to the number of stem cells and the number of progenitors in [3, electronic supplementary material, table S2] results in a poor 
agreement with the ratio of the number of stem cells to the number of progenitors. 



material, table S2]. For this purpose, we simultaneously fitted 
the analytical solution of (2.1) either for Ci and C1/C2 (the 
number of stem cells and the ratio of the number of stem 
cells to the number of progenitors), depicted in figure 6, or 
for Ci and C2 (the number of stem cells and the number 
of progenitors), shown in figure 7. Fitting was performed 
by using the NonlinearModelFit procedure of Mathematica 
V. 9 to numerically minimize the sum of squared residuals 
weighted by the inverse-square of the data points standard 
error of the mean as documented in the Mathematica reference 
for fitting data involving measurement errors. For numeri- 
cal minimization, the random search method was chosen, 
which resulted in the highest value among all available 
minimization methods. 

Fitting the solution for Ci{t) and Ci{t)/c2{t) of (2.1) to the 
stated data results in the values 

(laiOi -l)pi = -5.16 X 10-^h-\ ] 
KOi2(l-ai)pi = 171xlO-^h-\ > (5.1) 
d2 - (2^2 ^2 -1)P2 = 5.41 X 10"^ h-^ J 

Although this parameter set is in good agreement with the 
data of reference [3] for the age-related decay of the stem 
cell compartment and the dynamics on the ratio of the 
number of stem cells and the number of progenitors, the fit 



to the decay of the progenitor compartment is not good 
(figure 6). Conversely, fitting the solution for Ci{t) and 
C2{t) to the data of [3, electronic supplementary material, 
table S2], displays a good agreement of the stem cell and 
the progenitor compartment and a poor agreement of their 
ratio (figure 7). 

It appears that the decay of the stem cell compartment 
involves a saturation effect for late time points, which cannot 
be reproduced by our linear model. More precisely, the 
solution of (2.1) for the number of stem cells (ci) is an exponen- 
tially declining curve and fitting this curve to the data 
estimates five stem cells remaining at 2 years of age, whereas 
this number was measured to be 320. We hypothesize that 
this saturation is caused by either a feedback mechanism on 
stem cells that induces their quiescence with increasing age 
or by the existence of a mixture of two populations with one 
population performing adult neurogenesis and a quiescent 
one. Moreover, physiological conditions of Ui, pi, Oi and k, 
i.e. Ui ^ 0.5, 1/pi < 2 years, 0i > 0.1 and k > 0.5 contradict 
the restrictions imposed by (5.1), potentially because of the 
missing saturation effect that cannot be explained with our 
model. Thus, our current model indicates that there are some 
novel aspects in adult neurogenesis required to explain 
obtained experimental data such as this saturation effect. The 




Figure 8. Responses to an infinitesimal increase in stem cells fraction of self-renewal and proliferation rate, respectively, on the number of neurons, [a) Effect of an increased 
fraction of self-renewal with and without death of neurons, respectively, (b) Effect of an increased proliferation rate with and without death of neurons, respectively. 



proposed explanation of how saturation could be achieved 
should be the subject of future experimental validation. 



5.2. Simulations 

The mathematical analysis conducted in §4 depends on 
assumption 3.2. Because the purpose of this section is to 
extend our analysis to the effects of altered stem cell par- 
ameters on the number of neurons (q) and the number of 
BrdU-labelled cells {li{t,T)) for r > 0, we use a parameter set 
satisfying assumption 3.2 for our numerical investigations. 
Unless otherwise stated, we set ai = 0.55, = 0.7, = 1, 
K = 0.6, fl2 0.7, O2 0.4, p2 = 2.5, d2 = 0.1, ps 1.5, 0.4, 
i = 10, = 0.05, ds = 0.05, Ci(0) = 10 000, C2(0) = 5000, 
C3(0) = 15 000, C4(0) = 350 000, C5(0) = 100 000 to numeri- 
cally solve (2.1) by using the ND Solve framework of the 
Mathematica. Note that the stated parameter set does not 
include any time units, because a choice of physiological 
parameters was not feasible as stated in §5.1. Thus, there 
are no time units in any figure using this parameter set. 

At first, we investigate the effect of altered stem cell par- 
ameters on the number of existing neurons at time t after the 
change of stem cell dynamics, i.e. dpC4{t;p) for p G {ai, pi, 61]. 

An increase of Ui, which increases the proportion of 
symmetric stem cell divisions at the expense of asymme- 
tric divisions, displays the same qualitative progression on 
the number of neurons as the effect of an increase in Ui 
on the number of progenitors: initially, the neuron count is 
decreased, followed by an always ongoing period displaying 
an increase, but our simulations suggest that the magnitude 
of the increase weakens with time, i.e. \\m t^oo^ aiC^{t', cii) = 0. 
Interestingly, the existence of this weakening depends on the 
decay rate of neurons: if there is no decay {d/^ = 0), then we 
find that da^C4^{t;ai) converges to a positive value (figure 8a). 
Hence, the magnitude of the increase is not declining with 
time, if neurons are not allowed to decay. 

Increasing the proliferation rate pi of stem cells also affects 
the number of neurons in the same way as the number of 
progenitors, so that in the initial period, the neuron count 
is increased, in a subsequent period, the neuron count is 
decreased and the magnitude of the effect declines with 
time. Interestingly, if we assume, in this case, that neurons 
do not decay, the two-phase progression is lost and, the 
number of neurons is always increased, but the magnitude 
of the increase still weakens with time (figure 8h). 

Increasing the probability of stem cell divisions 61 dis- 
plays an increase in the number of neurons. If d^ > 0, our 
simulations indicate that the magnitude of the increase con- 
verges to zero and if d4 = 0 that it converges to a positive 
value (no figure is shown). 



Our analysis indicates that the effect of altered stem cell 
dynamics on the number of neurons depends on the death 
rate of neurons. Thus, investigating whether neurons have 
the ability to decay and if so determining their decay rate is 
vital in order to understand the short- and long-term impact 
of altered stem cell dynamics on the number of neurons in 
the dentate gyrus. 

Next, we turn our investigations to the quantity //(f*, r), 
the number of BrdU-labelled cells of cellular compartment i, 
where BrdU was applied at a particular time point t = and 
^ = 0 corresponds to the time point of the KO. The effect of an 
increase of stem cells fraction of self-renewal on BrdU- 
labelled cells depends on the time point f*. If BrdU was 
given directly after the fraction of self-renewal was increased, 
i.e. f* = 0, the number of BrdU-labelled progenitors ^(O/t) 
shows a two-phase progression with an initial decrease and 
subsequent increase. The same holds true for the number of 
BrdU-labelled neurons UiO, r) (figure 9, top row). As the 
time between the change of dynamics and the BrdU 
administration increases, the first phase that shows a decrease 
in the number of progenitors and neurons becomes shorter 
(figure 9, middle row) until there remains only a one-phase 
progression, with increased numbers of BrdU-labelled pro- 
genitors and neurons at every time point r after BrdU was 
given (figure 9, bottom row). 

Similar to an increased fraction of self-renewal, the effect of 
an increased stem cell proliferation rate depends on the time 
between the change of d3mamics and the BrdU administration. 
If f * = 0, then BrdU-labelled progenitors and neurons display 
a two-phase progression with increased numbers in the initial 
phase and decreased numbers subsequently (figure 10, top 
row). Interestingly, the two-phase progression after BrdU-label- 
ling gains an additional phase as increases, displaying now 
an initial decrease, a subsequent increase and again later on a 
decrease in the number of BrdU-labelled progenitors and 
neurons (figure 10, middle row). Increasing further, the 
three-phase progression is lost, and then the number of labelled 
progenitors and neurons is decreased at any time r after BrdU 
was given (figure 10, bottom row). 

6. Discussion and conclusion 

We have established a mathematical model of adult hippo- 
campal neurogenesis based on experimental data. Although 
we consider a basic model not accounting for any feedback 
mechanisms or a spatial component, we demonstrate that 
modifying the dynamics of adult neural stem cells, which 
corresponds to inducing a stem-cell-targeting KO, exhibits a 
rich variety of effects owing to the high complexity of the 
hippocampal neurogenic niche. 




Our investigation shows that observed differences in cell 
numbers owing to altered stem cell dynamics depend not 
only on the alteration that was induced by a particular KO, 
but also on the time at which cell counts were measured. 
Therefore, it is necessary to perform measurements at mul- 
tiple time points in order to draw conclusions from KO 
experiments. Moreover, we find that cell numbers of cells at 
different differentiation stages may respond qualitatively 
different to altered stem cell dynamics. Additionally, this 
response may also be different between the total number of 
cells and the number of BrdU-labelled cells of a given cellular 
compartment and also depends on the time point of BrdU 
administration after the KO. Thus, labelling cells with BrdU 
does not generate a subset of the neurogenesis system that 
is sufficient to perform an analysis of the impact of altered 
stem cell dynamics. The reason for this is that changing 
stem cell dynamics influences either the number of cell div- 
isions or the ratio of symmetric to asymmetric divisions, 
which, in turn, affects the initial distribution of cell types 
among all BrdU-labelled cells. Furthermore, our reasoning 
proves that under the assumptions of our model, the effect 
of altered stem cell dynamics declines as time passes and 
that this decline is a result of the depletion of the stem cell 
pool. Evaluating differences in cell counts at late time 



points after a KO can therefore be used to test the notion of 
a declining stem cell population. 

With the help of our study, the available data on KO 
experiments can be revisited in a new theoretical context: 
in reference [6], a KO of the gene Notchl is analysed at 
multiple time-points after the KO. The counts of stem cells, 
neural progenitors and mature neurons decreased at late 
time points, whereas they slightly decreased or remain 
unchanged at earlier time points [6] (figure 2). The authors 
conclude that Notchl KO results in a decreased self-renewal 
of stem cells. Moreover, they infer from neurosphere assays 
that neural progenitors also exhibit decreased self-renewal 
in the case of Notchl KO. Their reasoning and data are 
consistent with our theoretical findings: decreasing the self- 
renewal of stem cells (cf. §4.2.1) results in decreased counts 
of stem cells at any time after the KO, whereas the pro- 
genitor cell counts exhibit an initial growth followed by a 
subsequent decline. Additionally, decreased self -renewal of 
progenitors is observed. This simultaneous decrease can 
compensate the initial increase of progenitor numbers 
caused by the decreased self-renewal of stem cells such that 
the number of progenitors is also decreased at any time 
after the KO. In another study of Furutachi et al. [8], it is 
reported that p57 is responsible for the quiescence of adult 
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Figure 10. Effect of an infinitesimally increased stem cell proliferation rate on BrdU-labelled stem cells, progenitors and neurons at time r after BrdU-labelling, 
where BrdU was given at time f = 0 (a), t = 5.5 (b) and f = 10 (c). 



neural stem cells in the hippocampus: a KO of p57 resulted in 
a decrease of the number of BrdU-label-retaining stem cells at 
60 days after the KO (fig. 2k of [8]), whereas the progenitor 
count initially displayed an increase and a subsequent 
decrease (figs 4ib and 6/ of [8]). In the context of our model, 
the claim of the authors corresponds to the statement that 
the KO of p57 causes an increase of the stem cell proliferation 
rate pi. As stated in §4.2.2, an increase of pi results in an 
initial increase and subsequent decrease in the number of 
progenitors. Thus, the data presented and the authors' expla- 
nation that p57 KO decreases the quiescence of stem cells is 
consistent with the proposed mathematical model. 

Our modelling approach qualitatively describes the 
effects of altered stem cell d3mamics on cell counts and on 
BrdU-labelled cells. In order to calculate the magnitude of 
these effects and the time points at which positive and negative 
effects are to be expected, an estimation of model parameters 
accompanied by a possible data-driven modification of the 
proposed model has to be performed. Thus, a thorough charac- 
terization of wild-t3^e neurogenesis is needed in order to 
evaluate KO experiments within a mathematical framework 
and to draw conclusions about the underlying d3mamics of 
adult hippocampal neurogenesis. 
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Appendix A 

A.I. Proof of lemma 3.1 

It holds dci/dt < 0, if and only if aiOi < 1/2. To analyse 
the decay conditions of the progenitor compartment C2, we 
introduce the abbreviations 

q:= K6il{l-ai)pi, \ (Al) 
r:=d2-{2a2e2-l)p2. J 

First, let us note that p <0 and r > 0 is a necessary 
though not sufficient condition for 6.62/ dt < 0. Furthermore, 
it can be seen that the quantity C1/C2 satisfies a Riccati 
differential equation: 

dtWt))-c2it)V^' '^W' ^^^^ 

Because q > 0, the right-hand side of (A 2) describes a 
downward opening parabola with roots 0 and {p-\-r)/q. If 
p-\-r>0, the state {p-\-r)/q is globally asymptotically 
stable. Otherwise, 0 is globally asymptotically stable. In 



Table 2. Coefficients of dpci(t:p) for p G {a-i, pi, 6^}. 



E3 



fli a = (202^2 - 1)P2 - (2<Ji 01 - Dpi - (/2 a = -((/4 + (2(Ji 0i - 

T.-^K^P^n.-fli.) y=Mi^W-.M\-M+l-M 

Pi a = (202 ^2 - 1)P2 - (2fli 01 - l)pi - rf2 a = -((/4 + (2<Ji 0i - Dpi) 

/3 = K0i2(l - j8 = 0i(1 - k)2(1 - + 1 - 01 

y = K0i2(l - a,)(2aA - DPi 7 = HaA - 1)Pi(0i(1 - k)2(1 - fli) + 1 - 6,) 

a = (2fl2 02 -1) P2 - (2fli 0i - a = -((/4 T (2fli 0i - l) /?i^ 

/3 = k2(1 - fli)pi /8 = ((1 - k)2(1 - fli) - Dpi 

7 = Ke^M^ - a{|a^p^ T. =. 2fliX(0i(.l.-.K)2(1. .-.^ 

both cases, the steady state is approached monotonically. Fur- derivatives appearing in (A 9) have to be understood as par- 

thermore, p < 0 and q, r > 0 imply tial derivatives with respect to the considered parameter 

„ 1 ^ ^ p e jfli. Pi, f?i). It follows that 

< (^-3) 

1 1 / = (201^1 

and hence If y = C2, we have 

^ < L implies ^ < ^ for all t > 0. (A 4) S = '<ei2(l - a,)p^, 

'^2(0) 1 C2{t) q h = {7X1262 -\)p2-d2, 

Note from (2.1) that dc2/df < 0 holds if and only if 

and if 1/ = C4, it holds 

ci(f) ^ r 

^ = (ei(l-K)2(l-fli) + l-f?i)pi, 

Taken together with (A 4), we see that dci /df < 0 for all 
f > 0 if and only if explicit values of a, /3 and y appearing in (A 9) are sum- 

marized in table 2. 

^ Proof of lemma A.l. Define Xp{t) := dpX{t;p) and yp{t) := 

g dpy{t;p). Because x(0) and i/(0) are independent of p, we have 

Xp(0)=yp(0) = 0, (A 10) 

A.2. Derivatives with respea to parameters and it follows from symmetry of second partial derivatives 

Here, we derive the identities (4.4) -(4.5) and the parameter- that 

dependent values for a, p and y in the latter equation. For ^j^Xp , 

this purpose, we need lemma A.l, which we prove at the "dT^^^ — / (0^(0 

end of the section. j 

^(f) = g'(p)xit) + g{p)xpit) + hy,{t). 

Lemma A.l. Let p,h ^R, f, g real-valued, non-zero functions 

and consider the system Hence, 

f-nm. I ^^^^ ^(^)=/>-/'<f)-''<"5l=A.). 

^ = g{p)x{t) + hy{t), J The assumption x(0) 7^ 0 leads to {Xp/x){Q) = 0 by (A 10) and 

solving the above differential equation yields 
together with non-zero initial conditions x(0) and i/(0), which are ^ 
independent of p. The solution then satisfies ^^^^^ ~ f {p)tx{t). 

dpx{t; p) = f{p)tx{t\ (A 7) Furthermore, 

where yp(0) ^ 

a = h-f{p), I3 = g'(p) and y = f{p)g{p). (A 9) ^(0) 

imply 

In order to obtain (4.4) and (4.5), we apply lemma A.l to Vv^ ^ (e^^ - l){ap + y) - ayt 

(4.3) with X = Ci, y G {C2, C4} and p G [ai, pi, Oi]. Thus, the ^(0 



b 



with a, p and y as stated in (A 9), because a differential has the solution 

equation of the form 

x\t) = ax(t) + P + yt, x{t) = ^2 ■ 

x{0) = 0, 
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